On irregular prime power divisors of the 
Bernoulli numbers 



Bernd C. Kellner 

Abstract 

Let Bn {n = 0, 1, 2, . . .) denote the usual n-th Bernoulh number. Let I be a positive 
even integer where / = 12 or / > 16. It is weh known that the numerator of the reduced 
quotient \Bi/l\ is a product of powers of irregular primes. Let {p,l) be an irregular 
pair with Bi/l ^ Bi^p^i/{1 + p — 1) (mod/). We show that for every r > 1 the 
congruence Bm.^/rrir = (mod /) has a unique solution where = I (mod p — 1) 
and I < nir < {p — l)/~^. The sequence {mr)r>i defines a p-adic integer X{'p,i) 
which is a zero of a certain p-adic zeta function (^p / originally defined by T. Kubota 
and H. W. Leopoldt. We show some properties of these functions and give some 
applications. Subsequently we give several computations of the (truncated) p-adic 
expansion of X{p,i) for irregular pairs {p,l) with p below 1000. 
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1 Introduction 



The classical Bernoulli numbers B^ are defined by 

oo „ 

^^ = VS„^, \z\<2Tr. 
— 1 ^-^ n! 

n=0 

These numbers are rational where i?„ = for odd n > 1 and {—l)2~^^Bn > for even 
n > 0. The Bernoulli numbers are connected with the Riemann zeta function (^(s) on 
the positive real axis by Euler's formula 

C{n) = --^^Bn, neN,2\7i; (1.1) 
2 n\ 

the functional equation of C(s) leads to 

((l-n) = for n G N, n>2. (1.2) 

n 

Let, as usual, (p denote Euler's totient function. The Kummer congruences state 
for n, m,p,r G N, n, m even, p prime with p — 1 \ n that 

(1-/-^)— = (1-p™-^)— (mod/) (1.3) 
n m 

when n = m (mod ip{p'')); see [9, Thm. 5, p. 239]. 



In 1850 E. E. Kummer [ I I] introduced two classes of odd primes, later called regular 
and irregular. An odd prime p is called regular if p does not divide the class number of 
the cyclotomic field Q(^p) where fip is the set of p-th roots of unity; otherwise irregular. 
Kummer proved that Fermat's last theorem (FLT) is true if the exponent is a regular 
prime. Kummer also gave an equivalent definition of irregularity concerning Bernoulli 
numbers. We recall the usual definition from [9, pp. 233-234]. 

Definition 1.1 Let p be an odd prime. The pair (p, I) is called an irregular pair if p 
divides the numerator of Bi where I is even and 2 < I < p — 3. The index of irregularity 
of p is defined to be 

i{p) '■= is an irregular pair : I = 2,4, ... ,p — 3}. 

The prime p is called a regular prime if i{p) = 0, otherwise an irregular prime. 

We introduce the following notations for rational numbers. If q is rational then we 
use the representation q = N/D where {N,D) = 1 and D > We define denom(g) = 
D resp. numer(g) = N for the denominator resp. the numerator of q. The notation 
m I q where m is a positive integer means that m \ numer(g'); we shall also write 
q = Q (mod m) in this case. 

For now, let n be an even positive integer. An elementary property of the Bernoulli 
numbers, independently discovered by T. Clausen [C] and K. G. C. von Staudt [20] in 
1840, is the following. The structure of the denominator of Bn is given by 

Bn + — G Z which implies denom(S„) = p. (1-4) 

p—l\n^ p~l\n 

A further result, often associated with the name of J. C. Adams, see [9, Prop. 15.2.4, 
p. 238] and Section 8 below, is that Bn/n is a p-integer for all primes p with p — l\ n. 
Therefore 

W p^"^"'' (1.5) 

p— Ifn 

divides numer(i?„); since this factor is cancelled in the numerator of Bn/n, we shall 
can it the trivial factor of By \Bn\ > 2(n/27re)'", see [9, Eq. (8), p. 232], and 
Table A.l, one can easily show that the numerator of \Bn/n\ equals 1 for 2 < n < 10 
and n = 14. Otherwise this numerator is a product of powers of irregular primes; this 
is a consequence of the Kummer congruences. The determination of irregular primes 
resp. irregular pairs is still a difficult task, see [2]. One can easily show that infinitely 
many irregular primes exist; for a short proof see Carlitz [4]. Unfortunately, the more 
difficult question of whether infinitely many regular primes exist is still open. However, 
calculations in [2] show that about 60% of all primes less than 12 million are regular 
which agree with an expected distribution proposed by Siegel [18]. 

The values of Bn and Bn/n for n < 20 are given in Table A.l, irregular pairs with 
p < 1000 in Table A. 3. For brevity we write B(n) = Bn/n; these are called the divided 
Bernoulli numbers. Throughout this paper all indices concerning Bernoulli numbers 
are even and p denotes an odd prime. 
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2 Preliminaries 

The definition of irregular pairs can be extended to irregular prime powers as was first 
proposed by the author [IJ, Section 2.5]. 

Definition 2.1 Let p be an odd prime and n, I be positive integers. A pair {p^l) is 
called an irregular pair of order n if | B{1) where / is even and 2 < / < (f{p"'). Define 

$5.7 := {(p, : I B{1), p is an odd prime, 2<l< ipip'^), 2 | 1} 

as the set of irregular pairs of order n. For a prime p the index of irregular pairs of 
order n is defined by 

inip) :=#{(P,0: (p,Og<"}- 

Define the map 

A„ : - <\ (P, I) ^ {P, I mod (^(p")) 

where x mod y gives the least nonnegative residue x' with < x' < y for positive 
integers x and y. Let (p, In) S ^^^^ and (p, 1^) S be irregular pairs of order n resp. 
m where n > m > 1. We say that {p,ln) is related to (j),lm) if = (mod (fip^)) 
holds. Note that "related" is not a symmetric relation. 

Remark 2.2 This definition includes the older Definition 1.1 in the case n = 1. There- 
fore one has ii{p) = i{p)- Let (p, I) G with n>\. The Kummer congruences (1.3) 
imply that p" | B[l + !^(p{p"')) for all ly G No- Conversely, if | B{m) for some even 
integer m, then there exists an irregular pair {p,l) € where I = m (mod ^{p^)) 
with I < m holds. 

The map A„ is well defined by the properties mentioned above. Let (p, Z„) E 
and (j),lm) £ where n > m > 1 and {p,ln) is related to (j),lm)- By applying the 
maps A„_i, A„_2, . . . , A^ one derives a chain of related irregular pairs of descending 
order 

ip,ln) G ip,ln-l) G (p,Z„-2) G (P, ^m) G (2.1) 

where 

Definition 2.3 For (p,/) G ^'^^ n > 1 define 

A(p,o = p-" [b{1 + ^(p-)) - S(/)) (mod p) 
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with < ^[p^i) < P- When ^(pj) = we cah ^(p^i) singular. For an irregular prime p 
set 

\ 0, Ap = 

with 

i{p) 

Ap = nA(p,,„), (p,/.)en". 

Then A(p) = 1 if and only if all A(p^i^) are nonsingular. 

We need a generalized form of the Kummer congruences which allows us to obtain 
most of the later results, see Carlitz [3, Thm. 3, p. 425], especially Fresnel [7, Cor. 6, 
p. 319]. 

Theorem 2.4 (Carlitz) Let k,m,n,p,r,uj m even, p prime with p — 1 ] m, and 
Lv = k{p{p"'). Then 

^ ('''V-l)^(l-p™+'^^-^)5(m + z/cc^) =0 (modp""). (2.2) 

Here we need a special version without Euler factors which p-adically shows the 
periodic behavior of the divided Bernoulli numbers. 

Corollary 2.5 Let {p,l) G ^n^, n > 1. Let k,m,r,uj G N, r > 1, and oj = kip{p^). 
For m = I + j(p{p^) with j > we have 

Yl + '^uj) = (mod (p™-\p"("-i))). 

Proof. Since {p,l) G "S/^^ we know that | B{1 + j<f{p"')) for ah j > 0. Hence, we 
can reduce Congruence (2.2) to (mod p"(''~^)) when multiplying it by One easily 
sees that all Euler factors in the sum of (2.2) vanish (mod p"""!). □ 

Proposition 2.7 below shows how to find irregular pairs of higher order. Beginning 
from an irregular pair of order n, we can characterize related irregular pairs of order 
n + 1 if they exist. First we need a lemma. 

Lemma 2.6 Let n be a positive integer and p be an odd prime. Let {au)u>o be a 
sequence of p-integers £ Q for all v gNq. If one has 

ay - 2au+i + au+2 = (mod p") (2.3) 

then the sequence is equidistant (mod p"). For ao ^ ai (mod p) the elements ao up to 
otpn-i run through all residues (mod p"). Then exactly one element = (mod p") 
exists with < s < p"' where s = — ao(ai — cto)~^ (mod p"). 
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Proof. We rewrite Congruence (2.3) to 

au - Ou+i = a^+i - (modp"), G Nq. 

One easily observes by induction on v that all elements are equidistant (mod p"). 
Let S = ai — ao (mod p"), then we obtain 

= ao + 5v (mod p"). 

The case ao ^ ai (mod p) provides that ao + resp. runs through all residues 
(mod p") for < < p" , since 5 is invertible (mod p"). Then exactly one element 
as exists with = a^ = ao + (5s (mod p") and < s < p". This finally gives 
s = — ao(5~^ = — ao(ai — ao)""*^ (mod p"). □ 

Proposition 2.7 Xei (p, /) € ^n'^? n > 1. Let t/ie sequence {aj)j>o satisfy 

aj^p-^Bil + Mp^)) (modp). 

For A(p = ai — ao (mod p) where < A(p^/) < p there exist three cases: 

(1) If A^pj-j = and ao ^ (modp), then there are no related irregular pairs of 
order n + 1 and higher, 

(2) If A^^pj-^ = and ao = (mod p), then all {p,l + utp^p'^)) G ^"+1 related 
irregular pairs of order n + 1 for z/ = 0, . . . ,p — 1, 

(3) If A^pj-^ ^ 0, then exactly one related irregular pair of order n + 1 exists. One 
has {p, I + S{p{p"')) G ^n+i < s < p where s = — aoA^^^^ (mod p). 

Proof. Let j > 0. Using Corollary 2.5 with r = 2, lo = ^{p^), and m = l+j<f{p"') > 2, 
we get 

^ P J(-l)''p""i^(m + z^w) = (modp) 

which is 

aj — 2aj+i + aj^2 = (mod p). 

This satisfies the conditions of Lemma 2.6. We obtain three cases: 

Case (1): We have ao = ai (modp) and ao ^ (modp). One observes that 
aj ^ (mod p) resp. p"+^ | B{1 + j<f{p^)) for all j > 0. Therefore, there cannot exist 
related irregular pairs of order n + 1. Also there cannot exist related irregular pairs of 
order r > n + 1, otherwise we would get a contradiction to (2.1). 

Case (2): We have ao = ai (modp) and ao = (modp). This yields aj = 
(mod p) resp. p""*"^ | 1B{1 + jtp{p'^)) for all j > 0. Hence, p related irregular pairs of 
order n + 1 exist where (p, / + i'ip{p'^)) G ^^n+i for i/ = 0, . . . ,p — 1. 

Case (3): We have ao ^ ai (modp). Lemma 2.6 provides exactly one element 
as = (mod p) with the desired properties. Hence, (p, / + s(/5(p")) is the only related 
irregular pair of order n + 1. □ 
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Remark 2.8 Vandiver [19] describes the result of the previous proposition for the case 
n = 1 and only for the first irregular primes 37, 59, and 67. For these primes Pollaczek 
[16] has calculated the indices s of the now called irregular pair of order two, but case 
p = 67 with s = 2 is incorrect, see column S2 of Table A. 3. This error was already 
noticed by Johnson [10] who has also determined all irregular pairs {p, I') of order 
two with p below 8000. Wagstaff [21] has extended calculations of irregular pairs, 
indices s, and associated cyclotomic invariants up to p < 125 000. He also checked 
that FLT is true for all such exponents p in that range. Finally, Buhler, Crandall, 
Ernvall, Metsankyla, and Shokrollahi [2] have extended calculations of irregular pairs 
and associated cyclotomic invariants up to p < 12 000000. For all these irregular pairs 
{p, I) in that range ^{p,i) 7^ is always valid which ensures that each time there is only 
one related irregular pair {p,l') of order two. Hence i2{p) = i{p) for these irregular 
primes p. One has to notice that always {p,l) 7^ ip,l')- So far, no irregular pair {p,l) 
has been found with p'^ \ B{1). 

Using Proposition 2.7 one can successively find irregular pairs of higher order. We 
can easily extend the result starting from an irregular pair (p, I) G and requiring 
that I > n to obtain a related irregular pair (p, /') G ^fn- 

Proposition 2.9 Let {p,l) G ^n^, n > 1. Suppose that I > n. Let the sequence 
(aj)j>o satisfy 

aj = p-'^B{l + jV(p")) (mod p"). 
If Af^pj-^ ^ 0, then there is exactly one related irregular pair of order 2n 

ip,l + S^{pn)&^Tn 

with < s < where s = — ao(ai — ao)~^ (modp"). Correspondingly, there also 
exists exactly one related irregular pair of order n + u 

ip,i+sMpn)^K+u 

for each u = 1, . . . ,n — 1 with < Si, < p'^ where Sy = s (mod p^). 

Proof. Let j > 0. Using Corollary 2.5 again with r = 2, u; = y:'{p"'), and m = 
I + j(p{p'^) > n yields 

Oj — 2aj-|_i + aj+2 = (mod p"). 

If 7^ 0, then Lemma 2.6 provides exactly one element = (mod p") with 

< s < where s = —00(0^1 ~ oq)"^ (mod p"). Therefore, (p, I + s(/7(p'")) is the only 
related irregular pair of order 2n. Similarly, regarding the congruences above (mod p'^) 
instead of (mod p") for u = 1, . . . ,n — 1 yields the proposed properties. □ 

Finally, we start from an irregular pair (p, /) G where we have to suppose that 

1 > {r — l)n with r > 2 to obtain a related irregular pair (p, I') G ^'f^. 
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Proposition 2.10 Let {p,l) G ^n^, n > 1. Let r > 1 be an integer and suppose that 
I > (r — l)n. Let the sequence (aj)j>o satisfy 

aj = p-''B{l + jv?(p")) (mod p^'^"^)"). 

Then this sequence satisfies for all j > 0: 

The elements ao up to a^-i induce the entire sequence {aj)j>Q. Elements with = 
(mod where < s < p^^'~^^'^ provide related irregular pairs of order rn with 

{p,l + sip{p^)) G ^rn- ^(p,i) 7^ '^"■^ elements ao up to a^-i are equidistant 
(mod p^''"^^"), i/ien i/iere is exactly one related irregular pair of order rn 



^ rn 



{p,l + s^{p^))e% 

with < s < where s = — ao(Q^i ~ cko)~^ (mod p^^'"-^)'"). Correspondingly, there 

exists exactly one related irregular pair of order n + k 

{p^l + s,^{pn)^n\k 
for each k = 1, . . . , (r — l)n — 1 with < si^ < p^ where Sk = s (mod p^). 

Proof. Let j > 0. Clearly, by Corollary 2.5 with r > 1, u; = ifip""), and m = 
I + j'-p{p'^) > (r — l)n, we have 



^ ('■V-ira,+,- ^ (modp(^-i)"). 



y=0 

This induces the whole sequence (aj)j>o by 

r-l 

Otr+j 



g M (mod p(^-i)"). (2.4) 

Among all elements with < s < an element = (mod provides 

a related irregular pair (p, / + sy3(p")) € ^'J." of order rn. 

Now we assume that ^(p^i) 7^ and the first elements ao up to ar-i are equidistant 
(mod We show that this property transfers to the entire sequence. Let 7,5 

be integers. It is easily seen for r > 1 that 

jZ{'){-W{l + H=^- (2.5) 

Consider 7 = ao (mod p(^'~^)") and 5 = ai — ao (mod p(^~^)") where p \ 5 hy assump- 
tion. Combining (2.4) and (2.5) yields 

r-l . . 

^_iy+^Y,[){-W{ao + 5v) = ao + 5r = ar {mod p^'-^>) (2.6) 
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which shows inductively that all successive elements aj with j > r are equidistant 
(mod 

p{r-i)ny Since 6 is invertible (mod ^)"), exactly one solution exists with 
= as = ao+Ss (mod where s = — Qo(ai — ao)~^ (mod p^^"^)""). Using similar 

arguments, Congruences (2.4) and (2.6) are also valid (mod p^) for k = 1, . . . , (r — 
l)n — 1 which provides for each k a unique solution Sk with < Sk < where 
Sk = s (mod p^) holds. □ 

In [II, pp. 125-130] several examples and calculations are given which use the 
previous propositions. These results are reprinted in Table A. 4 and following. It will 
turn out later that calculations can be further simplified. Because of the rapid growth 
of indices, it is useful to write indices of irregular pairs of higher order p-adically. 

Definition 2.11 Let {p,l) G n > 1. We write 

n 

{p,si,S2,...,Sn) ^'^'n where I = Su (fip""^) 

for the p-adic notation of (p, /) with < Sj^ < p for v = 1, . . . , n and 2 < si < p — 3, 
2 I si. The corresponding set is denoted as '^n^, the map corresponding to A„ is given 

by 

A.„ : '^'^^i^'i'", {p,Si,S2,...,Sn,Sn+l) ^ {p,Si,S2,...,Sn). 

The pair (p, /) and the element {p, si,S2, ■ ■ ■ , Sn) are called associated. 

Remark 2.12 The definition of ^J^'^ means that we have = for n = 1. For 
n > 2 we can define a map "^"'^ —>■ '^n^, (p, I) ^ (p, si, . . . , s„) where the Sk are uniquely 
determined by the p-adic representation 

n-2 

l = Sl + {p-l)s, S = ^S„+2P'', 0<S^+2<P 

u=0 

and by si = / (mod p — 1) with 2 < si < p — 3. If Sfc = with k > 2 then there is an 
irregular pair (p, 1^) of order k with 

(p,/fc)e^f and (p,/fc)G^L-i- 

Note that (p, si, S2, ■ ■ ■ , Sn) is also called an irregular pair with (si, S2, . . . , s„) as the 
second parameter given p-adically. 



3 Main results 

Theorem 3.1 Let {p,li) € ^f''. If A^^pi^-^ ^ then for each n> 1 there exists exactly 
one related irregular pair of order n. There is a unique sequence (/n)n>i resp. {sn)n>i 
with 

(p. In) G 'I-r resp. (p, si , . . . , s„) G ^'^^^ 

and 

h < h < h ^ ■ ■ ■ , lim In = oo- 
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Moreover one has 

^{pM = ^{pM = \pM = ■■■ ■ 

IfA{p) = 1 then 

i{p) = i2{p) = h{p) = ... . 

Theorem 3.2 Let {p,ln) £ "^n^, 7i > 1 with A{p^i„) = 0. Then there are two cases: 

(1) (p, In) ^ ^n+i- There are no related irregular pairs of order n + 1 and higher, 

(2) {p,ln) G '^'n+i- There exist p related irregular pairs of order n + 1 where 
{p, /„+ij) G 1 with A(p^i^^^ ^) = and l^+ij = In+j'fip^) for j = 0, . . . ,p-l. 

The property of A(pj), whether ^{p^i) vanishes or not, is passed on to ah related 
irregular pairs of higher order. The case of a singular A(j, ;) would possibly imply a 
strange behavior without any regularity. By calculation in [2] up to p < 12 000 000, no 
such A{pj) was found. The following diagram illustrates both cases. 

Diagram 3.3 




Here a vertical line indicates that {p, In) G ^'J^'' n ^'JJ;!,.]^ happens. On the left side we 
then have p related irregular pairs of order n + 1 which are represented by branches. In 
this case the corresponding Bernoulli number B[ln) decides whether there exist further 
branches or they stop. Instead of n the order of the p-power must be at least n + 1. 
This also means that an associated irregular pair {p,si,. . . ,Sn+i) G "^n+i must have 
a zero s^+i = each time. In contrast to, the right side shows that in the case of 
^{p,h) 7^ there is only one related irregular pair of each higher order. If A{p) = 1 
then higher powers p'^ are equally distributed among the numerators of B{n). For each 
irregular pair considered, there exists exactly one index nk^u with n^^u = nQ^^ + kif{p'^), 
/c G No in the disjoint intervals 

ik^{pn,(.k+i)^{pn) 

for which p'^ \ B{nk^u) is valid. 

In [11, pp. 128-130] irregular pairs of order 10 were calculated for all irregular 
primes p < 1000. These results are reprinted in Table A. 3. In this table only one 
irregular pair has a zero in its p-adic notation: 

(157, 62, 40, 145, 67, 29, 69, 0, 87, 89, 21) G ^i'J. 
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Hence, one has with a relatively small index that 

(157,6 557686 520 486) G n ^'"7". 

It seems that these zeros can be viewed as exceptional; see also Table A. 2. It would 
be of interest to investigate in which regions such indices may occur. This could 
explain why no irregular pair (p, /) G ^^"^ D has yet been found, because these 
regions are beyond present calculations. Here we have index 12 000 000 in [2] against 
index 6 557686 520486. Because of the rare occurrence of zeros one can expect that 
{p, I) G n ^'Y I'esp- p"^ I ^(0 will not happen often. 

Another phenomenon is the occurrence of huge irregular prime factors. Wagstaff 
[22] has completely factored the numerators of the Bernoulli numbers Bn with index 
up to n = 152. Most of these irregular prime factors are large numbers, the greatest 
factors have 30 up to 100 digits. 

Finally, summarizing all facts together, the property A(p) = 1 can be viewed as a 
structural property of the Bernoulli numbers. This leads us to the following conjecture 
named by the author A-Conjecture. 

Conjecture 3.4 (A-Conjecture) For all irregular primes p the following properties, 
which are equivalent, hold: 

(1) ^[p^i) is nonsingular for all irregular pairs {p,l) G ^1", 

(2) A{p) = 1, 

(3) i{p) =i2{p) =i3{p) = ■■■■ 

Assuming the A-Conjecture one can also prove the existence of infinitely many 
irregular primes using only information about the numerators of divided Bernoulli 
numbers, see [11, Satz 2.8.2, p. 87]. Now we give the proofs of the theorems above. 

Proposition 3.5 Let {p,ln) G ^n^j n>l with A(p^;^-) 7^ 0. Then there is exactly one 
related irregular pair {p,ln+i) G "^"+1 ^(p,/n) = ^(p,/„+i)- 

Proof. We write A„ = A(^p i^y Note that /„ > 2 and p > 3. Define the sequence 
("j)i>o by 

aj=p~^B{ln+j^{pn) (modp2). 

Using Corollary 2.5 with r = 3, uj = ip{p"'), and m = In + j^{p^) > 2, we have for 
n(r - 1) > 2 that 

aj — 3aj+i + 3aj+2 — Oj+s = (mod p^). 

Taking differences with Pj = aj+i — aj yields 

(3j -2pj+i + Pj+2 = (modp2). 

By Lemma 2.6 the sequence {f3j)j>o is equidistant (mod p^). Proposition 2.7 shows that 
the sequence {cej)j>o is equidistant (mod p). By definition we have (3j = A„ (mod p). 
Therefore, we can choose suitable 7, (5 G Z so that 

Uj+i - Uj = /3j = A„ + p (7 + j6) (mod p^). 
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This yields 

i-i 

Oj = ap + y^(An + p (7 + 1^5)) 

From Proposition 2.7 we have 

s = -aoA~^ (modp), < s < p. (3.2) 

With In+i = In + sip{p"') we obtain the unique related irregular pair {p,ln+i) G ^n+i 
of order n + 1. As a consequence of Lemma 2.6, we observe that 

a^+jp = (modp). 

Thus we obtain a sequence ia'j)j>o defined by 

q;. ^ as+jp/p ^ + (mod p) 

which we can use to determine related irregular pairs of order n + 2 using Proposition 
2.7. By definition we have 

A„+i = a[ - Uq (mod p). 

It follows from (3.1) that 

pAn+i = p {a[ - a'p) = tts+p - as 

= pAn (modp^), 

since 

+ p\ fs\ 1 



2 J \2J 2^(^^ + 2.-1). 
Finally, we obtain the proposed equation A„+i = A„. □ 
Proof of Theorem 3.1. Using Proposition 3.5 with induction on n provides 

A(p,ii) = A(p_;2) = A(p,i3) = . . . 

with exactly one related irregular pair of order n 

{P, In) e resp. (p, si, . . . , Sn) G K'- 

The latter pair is given by Definition 2.11. Proposition 2.7 shows for each step n that 

In+l =ln + Sn V^(p"), < Sn < P (3.3) 
11 



with a suitable integer s^- This ensures that li < I2 < ^3 < ... as an increasing 
sequence ilj)j>i- Clearly, this sequence is not eventually constant, because | B{Iyi) 
with < \B{ln)\ < 00. Therefore lim„^oo = 00. 

Starting with (pji) G "^i^ this provides a unique sequence {lj)j>i- If we have 
another irregular pair {p,l'i) € with A(p ;/) 7^ and {p,li) 7^ iPjl'i) then 

/ /fc for all j, A; e N, 

because l'^ ^ Zi (mod and l[ = I'j (mod ip{p)) resp. Zi = (mod (p{p)) by (3.3). 

If A{p) = 1 then for each of the i{p) irregular pairs {p,h,u), v = 1, . . . ,i(p) there 
exists exactly one related irregular pair of higher order. Finally, it follows that i(p) = 
*2(p) = ^3(p) = • • • and so on. □ 

Proof of Theorem 3.2. Clearly, the (non-) existence of related irregular pairs in 
the case (1) resp. (2) is given by Proposition 2.7 case (1) resp. (2). Hence, we only have 
to show the remaining part of case (2). In this case we have p related irregular pairs 
(p, ln^\,v) = ip-, In + ^'■pip^)) S ^n+i Order n + 1 with = 0, . . . ,p — 1. Although 
A(p = we can use Proposition 3.5 by modifying essential steps. We then have 

aj = I3j = (mod p) . 

Congruence (3.2) must be replaced by 

s = 0,...,p-l (3.2') 

since A(p ^^-j = yields p values of s. It follows that 

^Wn) = = for z^ = 0,...,p-l. □ 



4 A p-adic view 

Let Zp be the ring of p-adic integers and Qp be the field of p-adic numbers. The 
ultrametric absolute value | |p is defined by \x\p = p-"^'^p^ on Qp. Define | |oo as the 
standard norm on Q. For n > 1 we define -i/'n : Zp — > Z giving the projection of Zp onto 
the set [0,^") n Z so that for x € Zp we have x — ipn{x) S p^Zp where < ipn{x) < p". 
We denote P as the set of the rational primes. Now, we shall use (1.2) to reformulate 
our results. Let (p, /) G "^"^ then 

_ c(i-0 - C(i-(^+P-i)) , , . 

A(p,o = (modp) 

with < A(p ;) < p. 

Theorem 4.1 Let (j>,l) £ with A(p 7^ 0. We define the sequence {ln)n>i recur- 
sively by li = I and, for n > 1, by 

Wl = ln + ^{p)i^n(^-^^^) = Z„+V^(p")Vlf%;^y 
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Then we have C{1 — 1^) G p'^'^p and consequently 

lim |C(1 — ln)\p = with In — > oo. 

n— >oo 

Proof. We rewrite our results using <^(1 — /„) = —B{ln). Theorem 3.1 provides one, 
and only one, sequence {ln)n>i with = I and (p. In) G '^n^ with /„ oo. This implies 
the p-adic convergence C,{\ — In) 0. Additionally, from Proposition 2.7 we have 

In+l = In + Sip{p"-) 

for each step where 

s = A(p^;) (mod p), 0<s<p. 

Rewriting the last congruence yields 

The rest follows by tpniap""^) = ^i(o) for a (i'Lp. □ 

These results can be also applied to the so-called p-adic zeta functions which were 
originally defined by T. Kubota and H. W. Leopoldt [1.3] in 1964; for a detailed theory 
see Koblitz [12, Chapter II]. We introduce some definitions. 

Definition 4.2 Let p be a prime. For n > 1 define 

Cp(l-n) := (l-p"-i)C(l-n) = -{I - f^~^) B{n). 

Define the p-adic zeta function for p > 5 and a fixed si G {2, 4, ... ,p — 3} by 

Cp,si : ^p, Cp,si(s) := Inn Cp(l - (si + (p - 1)*^.)) 

resp. for p > 2 and si = by 

Cp,o : ^p \ {0} ^ Qp, Cp,o(s) := lim Cp(l - {p - l)tu) 

for p-adic integers s by taking any sequence {tu)u>i of nonnegative integers resp. pos- 
itive integers in the case si = which p-adically converges to s. 

Remark 4.3 Let p > 2 and si > 0. The p-adic zeta function Cp,si{s) interpolates the 
zeta function Cp(l — n) at nonnegative integer values s by 

Cp,si('5) = Cp(l - n) 

where n = si (mod p — 1) and n = si + [p — \)s. 

Let p > 5 and si 7^ 0. The Kummer congruences (1.3) state for r > 1 that 

Cp, (s) = Cp, si (s') (mod p") 

when s = s' (mod p*""^) for nonnegative integers s and s' . Since Z is dense in Zp, the 
p-adic zeta function C,p^ extends uniquely, by means of the interpolation property and 
the Kummer congruences, to a continuous function on Zp; see [12, Thm. 8, p. 46]. 
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Definition 4.4 Let {p, I) € with ^(p^i) 7^ 0. Define a characteristic p-adic integer 
which contains all information about irregular pairs of higher order by 

where {s^)^^i is the sequence defined by Theorem 3.1 and I = si. 

Lemma 4.5 Let (p, I) G with A(p /) ^ 0. Let r = ordp B{1) and {p,si, . . . , Sr+i) S 
^r+i ^6 related irregular pair of order r + 1. Then 

•Sr+i ^{p,i) = -p~'' B{1) (mod p) 

with si = I, = for v = 2, . . . ,r, s^+i 7^ 0. If r = 1 then X{p,i) S otherwise 
X{pj) e ^'■"^Zp /or r > 2. 

Proof. Since r = ordp ^ 1) we have (p, /) € for all = l,...,r. Then 
Proposition 2.7 and Theorem 3.1 provide 

s^-p-'B{l)A~^]i^ (modp), {p,l + s^{pn)e^';U 

with < s < p. We have s ^ since ordp {p~'^B{l)) = 0. Let (p, si, . . . , s^+i) G ^r+i 
be the related irregular pair of order r + 1. Then we see that si = /, s^^i = s ^ 0, and 
3,^ = for z/ = 2, . . . , r. By Definition 4.4 we have 

= S2 + S3P + . . . + Sr+ip"""-^ + . . . . 

Case r = 1 yields S2 7^ and X(p, «) S ^p; otherwise case r > 2 implies that ordp X{p,i) = 
r - 1 and X(p,/) ^ p'"-^Zp. □ 

Theorem 4.6 Let {p,l) G ^f'^ twii/i A(p, /) 7^ 0. The p-adic zeta function Cp,iis) ho-s a 
unique zero at s = X(p.i)- 

Proof. Prom Theorem 3.1 we have a sequence {ln)n>i with I = li. In view of 
Definition 4.2, Theorem 4.1 also states that 

lim |C(1 — ^n)|p = with — > cxd. 

n— >oo 

We can transfer this result to the p-adic zeta function (^p j by the interpolation property. 
We see that p-adically 

lim 1^ = 1 + {p- l)x(p,i)- 

Since the function C,p^i is continuous, the p-adic integer X{p,i) is a zero of ("p,/- We shall 
show that this zero is unique. Assume that Cp.iiO = with some ^ E Zp. We can use 
the arguments given above in the opposite direction. Since C,p^i is continuous, there 
exists a sequence {l'n)n>i of positive integers with 

lim l'^ = l + {p- 1)^ and lim |Cp(l - /„)|p = 0. 
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We can choose a subsequence {ln)n>i of (Z^)n>i such that — l'^) G p^'^p- By use of 
the Kummer congruences, we construct the sequence {ln)n>i where In = I'n (mod ^{p^)) 
with I < In < ^{p"^)- Now we have (pjn) € ^n"^ for all n > 1. Since I = h = h, 
Theorem 3.1 shows that (/n)n>i = {ln)n>i- This implies that ^ = X(p,0- '-' 

Remark 4.7 Let (p, G 'ff'' with A(p^;) 7^ 0. The related irregular pair (p, si, . . . , Sr) 
€ ^'^'^'^ is a p-adic approximation of the zero X(p,i) of the p-adic zeta function (pj. For 
the first irregular primes 37, 59, and 67 elements of ^foo were calculated in [11, pp. 127- 
128]. These results are reprinted in Table A. 2. By Lemma 4.5 the statement f Bi is 
equivalent to the fact that X{p,i) is a unit in Zp. 

From now on, we shall assume the A-Conjecture. We shall see that the zeros X(p, i) 
play an important role in the representation of the Riemann zeta function at odd 
negative integer arguments. 

Theorem 4.8 Let n be an even positive integer. Under the assumption of the A- 
Conjecture we have 

'n\p -p-r p 



c(i-n)=(-i)t n ^ n 



l = n (mod p—1) 

Proof. Since both products above have only positive terms, the sign follows by (1.2). 
The first product describes the denominator of C(l — n) which is a consequence of (1.4) 
and (1.5). We have to show that the second product describes the unsigned numerator 
of (^(1 — n) which only consists of powers of irregular primes. Let p be an irregular 
prime divisor of C(l — n). From Remark 2.2 we have 

ordp i3(n)=r =^ n = Ir {mod (f{p^)) with (p, /r)€^^". 

The irregular pair {p,lr) of order r is related to some irregular pair {p,l) € with 
A(p,/) 7^ where I = Ir = n (mod ^p{p))- We also have by Definition 4.4 that 

n = lr = l + {p-l) X{p, I) (mod (p - l)/~^Zp) 
and equally by reduction that 

= X{p,l) (mod p''"^Zp). 

The last congruence is not valid (mod p'^l^p) by construction. Therefore 

n — I 



X(p,i) 



p—1 



p-{r-l) 



which provides, with an additional factor p, the desired p-power in the second product. 
Considering all irregular primes p which can appear, the second product equals the 
numerator of C(l — n) without sign. □ 

With some technical definitions we can combine both products of the theorem 
above. This yields a more accessible representation of the Riemann zeta function by 
means of p-adic analysis. 
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Theorem 4.9 Define = ^1" U (P x {0}) and set X(p,o) = for all p £F. Define 
p{l) = 1 — 2sign(/) = ±1 for I > 0. Let n be an even positive integer. Under the 
assumption of the A- Conjecture we have 



c(i-n)=(-i)t n 



\X{p,l) \ 



(p,i)e*o 

l = rL ( mo cl p — 1) 



Proof. We only have to consider case / = 0. Then we have p—1 \ n, p{0) = 1, and 
IX(p.o) ~ ^^\p ~ b^lp- '^^^ other case / > is already covered by Theorem 4.8. □ 

We shall give an interpretation of this formula above in Remark 4.17 by means 
of p-adic zeta functions. This generalization shows the significance to prove the A- 
Conjecture at all. 

Theorem 4.10 Let {p,l) G ^'i" with A(p^;) ^ 0. Let s,t £ Zp. Then a strong version 
of the Kummer congruences holds that 

\CpA^) -Cv,iit)\p = \pi.s-t)\p. 

Moreover 

-\p,l) (modpZp) for s^t 



Cp,lis) - Cp,l(t) 



p{s-t) 
and 

Cp,iis) = -P^{p,i) (mod p'^Zp). 

Thus, A(p^;) is closely associated with thep-adic zeta function Cp,/ in the nonsingular 
case A(p ;) ^ 0. We will prove this theorem later. 

Corollary 4.11 Let {p,l) S ^'x" with A(p ^ 0. The p-adic zeta function Cp,/(s) has 
a simple zero at s = X{p,i)- Moreover, for s £ Zp, 

Cp,i{s) =p{s- X{p,i))Cp,i{s) 

where ^(s) is a continuous function on Zp with ^{s) = — A^p (mod p7L^. Con- 
sequently 

\C,p,i{s)\p = |p(s -X(p,^))lp• 
PROOF. Theorem 4.6 shows that Qp^iis) has a unique zero at s = X(p,/)- We can use 
Theorem 4.10 to define 



CpM 



s + Xip,i), 



P{s-X{p,i)) 
CpA^) 



p s = X{p,l)- 
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Theorem 4.10 implies that Cp li^) = ~^{p,i) (mod pZp) for all s G Zp. Hence 
C* lis) has no zeros and consequently Cp,z(^) ^ simple zero at s = X(p,i)- Since 
Cp,i{s) is continuous on Zp and Cp exists at s = X{p,i)i t^ie function Cp 
is also continuous on Zp. Finally we obtain Cp,/(s) = p{s — X{p,i)) Cp and 
\Cp,i{s)\p = \pis - X{p,i))\p lCp,/(s)lp = \pis- X{p,i))\p- □ 

Definition 4.12 Let p be a prime. Define 

C^o : ^P - ^P. Cp,o{s) := - lim [1 - p^^iP-^)-^)E^MEzll 

for p-adic integers s by taking any sequence {tu)u>i of nonnegative integers which 
p-adically converges to s. 

Proposition 4.13 Let p be a prime. The function q(s) is continuous on Zp \ {0} 
and satisfies the following properties: 

C*o(o) = -i, C,o(s) =PsCp,o(s), seZp\{0} 

and 

lim \ps CpAs) - Cp,o{0)\p < 1- 

s— >0 

Moreover Cp, o(^) = (mod pZp) for s G Zp in case p > 2 and s G pZp in case 
p = 2. Additionally, if p = 2 then Cp, o('^) ~ ^ for s = 1 (mod pZp). 

Proof. From Definition 4.2 and Definition 4.12 we have for s G Zp \ {0} and any 
sequence {ti,)i,>i of positive integers which p-adically converges to s that 

psCpAs) = - iim(i-p*-(p-i)-i)p. f;-(^";) 

= -lim(l-p*^(^-i)-i)^^^ = C;,o(^). 

Moreover we have 

C*o(0) = -(l-p-i)pi?o/(p-l) = -l. 

By Clausen-von Staudt (1.4) we obtain Cpo('^) = ~^ (mod pZp) for s G Zp in case 
p > 2 and s G pZp in case p = 2. Additionally we have C2 o(^) ~ ^ ^'^'^ s G 1 + 2Z2, 
since Q oi^) = ~(1 ~ 2°) 2i?i = and i?„ = for all odd integers n > 1. We use the 
fact that Cp,o : Zp ^ Qp is a continuous function on Zp \ {0}; see [12, Thm. 8, p. 46]. 
Hence q(s) = ps Cp,o(s) is also continuous on Zp \ {0}. It remains to show that 

lim |psCp,o(s) -C;,o(0)lp < 1- 

s— >0 

We can choose a zero sequence (sjy)i/>i where its elements are arbitrary close to 0, say 
< \su\p < P~^ with some positive integer r. Then we deduce that 

psu Cp,o(si.) - Cp,o(0) = Cp,o{su) + 1 = (mod pZp) 

for all v > 1. This implies, together with the continuity ps (^p^q{s), the estimate 
above. □ 
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Corollary 4.14 Let p be a prime. The p-adic zeta function Cp,o{s) ^ols a simple pole 
at s = 0. Moreover, for s ^ Zp\ {0}, 

Cp,o{s) = Cp,o(.s)/{ps). 

Consequently 

\Cp,o{s)\p = \ps\~^ 
for s & Zp \ {0} in case p > 2 and s £ pZp \ {0} in case p = 2. 

Proof. By Proposition 4.13 we can write Cp,o{s) = Cp^Q{s)/{ps) for s e Zp \ {0}. 
Moreover linis^ops Cp,o(s) = C with |^ — Cp, o(^)lp < ^- This implies that the limit 
lim5_>o Cp, o(s) does not exist and that Cp,o(s) has a simple pole at s = 0. We have 
\Cp o('^)lp = 1 for s € \ {0} in case p > 2 and s G pZp \ {0} in case p = 2. In these 
cases we then obtain |Cp,o('5)|p = b'Slp"^- D 

Remark 4.15 One can even show further that Cp o(^) continuous on Zp where 
lim^^o Cp,o(^) ~ Cp.o(O); moreover Cp^oi^) satisfies the Kummer congruences. This can 
be derived, e.g., by means of p-adic integration, see Koblitz [12, pp. 42-46], or by using 
certain congruences of Carlitz [5]. 

Theorem 4.16 Define = ^-i" U (P x {0}) and set X{p,o) = for all peF. Define 
Sp,i{n) = {n — l)/{p — 1). Let n he an even positive integer. Then 

ic(i-n)u= n \^p.iM^w= n \cp,iM^w- 

pep {p,Oe*o 

l = n (mod p — 1) Z = n (mod p— 1) 

Under the assumption of the /^-Conjecture we have 

ic(i-?^)ioo= n b(«p,'W - x(p,/))ipy n b(«p-oH-x(p,o))ip^- 

(p,OG*\" p-l\n 

l = Ti {mod p — 1) 

Proof. Since n is an even positive integer, the product formula states that 

n ic(i-n)ip=i. 

pePu {cxd} 

By Definition 4.2 we have |C(1 — ?^)|p = \Cp,i{s)\p where / = n (mod p — 1) with 
< / < p — 1 and s = Spj{n) = {n — l)/{p — 1). Thus 

m-n)\'J=\{m-n)\p= n \CpAspAn))\p. 

peP pep 

l = n (mod p — 1) 

We have \Cp,i{s)\p = 1 for s G Zp if (p, I) ^ and / / 0. Hence 

n \cpA^p,i{^))\p = n \^pa^pA'^))\p ■ 

p6P (p,06*o 
! = n (mod p-1) l = n (mod p-1) 
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Next we assume the A-Conjecture. By Corollary 4.11 we then have for {p,l) G 
that \Cp,i{s)\p = \pis — X{p,i))\p foi' s € Zp. Without any assumption, Corollary 4.14 
shows that |Cp, oCs)!^ = = \p{s — X(p,o))\p^ for s € Zp \ {0} in case p > 2 and 

s G pZp \ {0} in case p = 2. Since n is even and S2,oin) = n, we finally obtain 

n \Cp,lisp,lin))\p = Y[ \p{spA{n) - X(p,i))\p n \Pi^P,oin') -X{p,o))\p^- 

(p,Oe*o (p.Oe*']^"^ p~l\n 

l = n (mod p-1) (^^j 

□ 



Remark 4.17 Assuming the A-Conjecture, the numerator of |C(1 — n)|oo with even 
n > is essentially described by simple zeros X{p,i) of p-adic zeta functions Cp, K*) 
where I = n (mod p — 1) and (p, /) G ^'^f'^. More precisely, by the variable substitution 
s = {n — I) / {p — 1), the term 

\p{s - X{p,i))\p^ 

is equal to p''' for some suitable r > which is the p power factor in the numerator of 
|C(1 — ri)\oo. The term above is induced by 

\Cp,i{s)\p = \p{s - X{p,i))\p =P~^\{s - X{p,i))\p- 

Since X(p,i) is a simple zero of C,p,i{s)^ \C,p,i{s)\p is mainly determined by a linear factor 
which one can also interpret as a distance between s and X{p,i)- 

Similar arguments can be applied to the denominator of |C(1 — n)|oo- Without any 
assumption, the denominator of |C(1 — n)|oo is essentially described by simple poles 
X(p,o) of p-adic zeta functions Cp,o(s) where p — l\n. Again, the term 

b(s -X(p,o))lp^ 

is equal to p"^ for some suitable r > which is the p power factor in the denominator 
of |C(1 — ?i-)|oo- This term is induced by 

ICp,0(s)|p= |p(s-X(p,0))lp^ = P\{S - X{p,Q))\p^ 

where Cp,o('5) has a simple pole at s = X(p, o)- 

Now we shall make some preparations to give later a proof of Theorem 4.10. 
Definition 4.18 Define the linear finite-difference operator V and its powers by 

for r > and any function / : Zp — > Zp. The series 

with coefficients G Zp where \ay\p is called a Mahler series which defines a 
continuous function / : TL^^TLp. 
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The following theorem of Mahler shows that the converse also holds, see [17, Thm. 1, 
p. 173]. Note that the sign (— 1)"^ depends on the definition of D. 

Theorem 4.19 (Mahler) Let f : ^ "Zp be a continuous function. Then f has a 
Mahler series 

where the coefficients are given by = {—\YT)^f{^) and \a^\p holds. 

Definition 4.20 Let {p,l) £ ^''f with A(p ;) 7^ 0. For n > 1 define the p-adic zeta 
function of order n by 

Cp,l,n{s) =P~'^Cp,l {i^n-l{X{pA)) +P'^~'^s), S £ Zp 

which for s = corresponds to the related irregular pair of order n. 
Proposition 4.21 Let {p,l) G with ^(pj) 7^ 0. For positive integers n, r we have 
V"^\p,iA^) = ^ (modj>"%), s^Zp (4.1) 

and 

Cp,Ln{s) = Cp,l,n{t) (modp'^Zp), s,tGZp (4.2) 
when s = t (mod p^Zp). 

Proof. First assume that s G Nq. In analogy to Corollary 2.5, we have to modify 
Theorem 2.4 in a similar way, since Cp, i,n(0) corresponds to the related irregular pair 
of order n. Rewriting (2.2) in the case k = 1 gives (4.1). Since Z is dense in Zp, we 
can extend (4.1) to values s G Zp by means of the interpolation property of Cp,/ resp. 
Cp,l,n- By the same arguments the Kummer congruences, given in Remark 4.3, can be 
extended to values in Zp. Let s = t (mod p^Zp) and write s' = tpn-iix^pj)) + p^~^s 
and t' = i;n-i{x{pj)) + P^'^t- Then s' = t' (mod p^+^^^Zp) and therefore (?,«(«') = 
Cp,K*') (mod p'"+''Zp). By Definition 4.20 this gives (4.2). ' □ 

Proposition 4.22 Let {p,l) G with ^{pj) 7^ 0. For 7i > 1 we have 

Cp,l,nis) = ^(p,l){sn+i- s) (modpZp), S G Zp 
where Sn+i is defined by X{p,i) = S2 + ssP + • • • • There exists the Mahler expansion 

Cp,l,n{^) = Cp,/,n(0) + S G Zp 

u>l 

with G Zp where Zy = (— l)'^p~"('^~i) X''^Cp,i,n(0) and z\ = —A(^pj^ (mod pZp). 
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Proof. Proposition 2.7 also works with aj = Cp,i,n{j) (mod pZp) for j € Nq. Since 
Cp, i,n(0) corresponds to the related irregular pair of order n, we then have 

Cp,i,n(0) = A(p ;) Sn+1 (mod pZp) 



where is given by Definition 4.4. A further consequence of Proposition 2.7, ex- 

^^Cp,/,n(s) = A(p^;) (modpZp) (4.3) 



tended to s € Zp, is that 



and furthermore 

Cp,i,n{s) = A(p^;) {sn+1 - s) (mod pZp). 

By construction Cp,/,n is continuous; this is also a consequence of (4.2). Theorem 
4.19 shows that Cp,i,n has a Mahler series with the coefficients = (— 1)"^ P'^Cp.i,n(0) 
for u > 0. Using (4.1) we see that G p'^^'^~'^^'Zp for v > 1. Thus, we can set 
for v > 1 to obtain the proposed series expansion above. From (4.3) 
we deduce that zi = — A^^^^^ (mod pZp). □ 

Corollary 4.23 Let {p,l) € with A(p^i) / 0. Let 

U>1 



6e i/ie Mahler expansion of (^pj^i given by Proposition 4-. 22. We have for r > 1 that 

r-l 



Cp,K«) = Cp,KO) + j (mod/Zp), s G 



Special cases are given by 



CpAo) = -Y.p'^zJ 



X{p,i) 

V 



and for r > 1 ; 



Cp,KO)^-j;p'^^.(^^^''^) (modp%) 



Proof. We rewrite the Mahler expansion above using C,p.i,i{s) = p ^C,p.i{s) which 
yields 



Cp,i{s) = CpM + ^-{'X '^'^p- 

U>1 



By the assumption A(p 7^ we have the zero X{p,i) of Cp,i- The finite sums follow 
easily (mod p'^Zp), since the coefficients p'^Zi, vanish for v >r. □ 

Hence, we can use Corollary 4.23 to verify calculations of the coefficients and of 
the zero X{p,i)- Now, we are almost ready to prove Theorem 4.10; we first recall the 
following lemma from ['", p. 227]. 
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Lemma 4.24 For k > 1 and p> <k < p^"*"^, we have 



Proof of Theorem 4.10. Let s,t eZp with s ^ t and set r = ordp(s — t)>0. Then 
we have by Kummer congruences that 



s = t (mod p^l^p) 
We show that 

Cp,i{s) - Cp,i{t) 



CpAs) = CpAi) (mod/+%). 



-A(p;) (mod pZp) for s^t 



(4.4) 



(4.5) 



Pis-t) 

which also imphes the converse to (4.4). Set n = min {ordp Cp./(*)) ordp(^pj{t)} where 
n > 1. By Definition 4.20 we rewrite s = V'n-i(X(p,o) +P^~^ s' and t = tljn-i{x{p,i)) + 
t' with some s',t' £ Zp. Set u = r + l — n>0. Note that ordp(s' — t') = u. By 
the imphcation in (4.4) we then have 

p""(Cp,/(s) - QpAt)) ^ Cp,iAs') - Cp,iAt') = (mod p"+iZp) 

with some integer c. On the other side, we can use the Mahler expansion of Cp./,n given 
by Proposition 4.22. We obtain 



Cp,l,n{s') - Cp,l,n{i') = Zi{s' -t')+p''z2 



+ p'''z. 



+ 



(mod p"+^Zp 



Since ordp(s' — t') = u, Lemma 4.24 shows that 



ord,; 



>u - [log k\ 



where logp is the real valued logarithm with base p. Since we have p > 5 we obtain the 
estimate 



ordp p 



n(fc-l) 



Zk 



>u + n{k-l)- [log k\>u + l 



for all k > 2. Hence, all terms vanish (mod p"+^) except for k = 1: 
Cp,i,nW) - Cp,iAt') = ^1 - i') (mod p"+iZp). 

By Proposition 4.22 and ordp(s' — t') = u we finally get 

zi (s' - t') = - A(p, {s' - t') = p^'c (mod p"+iZp) 

which shows that c ^ (mod p). Since p~^^^ {s — t) = s' — t', this deduces (4.5). 

Now, we have to determine the derivative of Cp.i- Taking any sequence {tiy)^>i with 
ty ^ s for all 1/ > 1 and limjy_^oo = s, (4.5) yields the derivative with 

Cp,iis) = -pA(j,^i) (mod/Zp). ^ 
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Remark 4.25 In general the converse of the Kummer congruences does not hold. The 
first nontrivial counterexample is given by p = 13 and Biq/16 — -B4/4 = —7 • 13^/2720. 
The only case where the divided Bernoulli numbers are equal is Bn/IA — B2/2 = 0. 

At the end of this section we revisit the p-adic zeta functions Cp, si in the case 
si 7^ as introduced by Definition 4.2. We can transfer some results to these functions. 
Moreover, we can state a formula equivalent to the Kummer congruences, but which 
involves values of p-adic zeta functions at smallest possible argument values. 

Proposition 4.26 Let p he a prime with p > 5 and si S {2,4, . . . ,p — 3} be fixed. 
Then 

VCp,sAs)=0 (modp%), sGZp. (4.6) 
There exists the Mahler expansion 

with Zy G Zp where Zy = (-!)>-'' r'''Cp,si(0). 

Proof. In analogy to Propositions 4.21 and 4.22 in the case n = 1, Congruence (4.6) 
is a consequence of Theorem 2.4. Since Cp, si is ^ continuous function on Zp, we obtain 
a Mahler expansion where the proposed coefficients follow by (4.6). □ 

Theorem 4.27 Let p he a prime with p > 5 and si G {2, 4, ... ,p — 3} be fixed. For 
r > we have 



Cp,si(s) = '^Cp,si{k) Tr^k{s) (mod/+^Zp), s G Zp 
fc=o 

where the polynomials Tj.^k : Zp ^ Zp with T^^k £ o-nd degT^. ^ = r are given by 

j—k 

Proof. We rewrite the Mahler expansion of Cp, si given in (4.7). Let s G Zp. For r > 
we obtain the finite expansion 

Cp,.i(s)=Cp,.i(0) + X]p^^.(') ^E(-ir^^Cp,.i(0)(') (modp'^+%). 

By Definition 4.18 we have 

^^Cp,.i(0) = ^ p^)(-l)%,.,(A:) =^ (^^)(-l)%,.,(A:). 
fc=0 ^ ^ fc=o ^ ^ 
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Reordering the finite sums and terms yields the definition of T^^k above. Since x \—^ 
which is a polynomial of degree j defines a function on Zp, the polynomials T^^k are 
functions on Zp. The coefficients of T^^k are rational, e.g., rewrite Tr^r{x) = (^) as a 
polynomial in x. We deduce that degT^^^ = r, because the term (^) which occurs only 
once gives the maximal degree. □ 

Corollary 4.28 Let n he an even positive integer and p he a prime where p — \ \ n. 
Define the integer si hy si= n (mod p—l) with < si < p—\. Set s = {n — si)/{p—\). 
For r > 1 we have 

r-l 

(1 - p^~')Bin) ^ ^^(1 - ^^+'=(^-1)^1)5(51 + Hp - 1)) Tr.i,k{s) (mod /Zp) 
with Tr^k O'S defined ahove. 

Proof. This is a reformulation of Theorem 4.27 using Definition 4.2 where si + s{p — 
l)=n. □ 



Remark 4.29 The case r = 1 of Corollary 4.28 reduces to a special case of the Kum- 
mer congruences (1.3) for r = 1. We explicitly give the cases r = 2,3, 4 of Corollary 
4.28. Note that si > 2 and p>5. Case r = 2: 

(1 - p"-^)S(n) = -{s - 1) (1 - p''~^)B{si) + s B{si +P-1) (mod /Zp). 

Case r = 3: 

(l-p"-i)B(n) ^ +l(s^-3s + 2){l-f^~')Bisi) 
-s{s-2)B{si+p-l) 

+ ^s{s-l)B{si + 2{p-l)) (modp%). 

Case r = 4: 

(1 -p'^-i)5(n) = -l(s3-6s2 + lls-6)(l-/^-i)^(si) 

6 

+ ^s(s2 -5s + 6)5(si +P-1) 

-is(s2-4s + 3)S(si + 2(p-l)) 

+ ls{s - l){s - 2) B{si + 3{p - 1)) (mod p%). 
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5 Algorithms 



Here we will give some algorithms for calculating irregular pairs of higher order as- 
suming we are in the nonsingular case A^^ 7^ 0. As a result of Theorem 3.1 and 
Theorem 4.1, one has first to calculate ^(pj), then the irregular pair of order n resp. 
the corresponding divided Bernoulli number provides the next related irregular pair of 
order n + 1. This is not practicable for higher orders say n > 3. Proposition 2.10 shows 
a more effective way to determine related irregular pairs of higher order. Starting from 
an irregular pair {p,l) € with n > 1 and requiring that I > [r — l)n with some 
r > 2 we can obtain a related irregular pair (p, V) € ^^n- If the corresponding sequence 
(aj)j>o is equidistant (mod then one can easily apply this proposition. If not, 

one has to calculate successively all elements Uj in order to find = (mod 
where < s < p^**"^)" which exists uniquely by the assumption ^(p.i) 7^ 0- To shorten 
these calculations, this search can be accomplished step by step, moving each time 
from a sequence (aj,fc)j>o to a sequence (aj,fc+i)j>o which are assigned to the irregular 
pair of order k resp. k + 1. 

Proposition 5.1 Let {p,l) G ^"'^ , n>l with A(p ;) 7^ 0. Let r, u he positive integers 
with r > 1 and u = {r — l)n. Assume that I > u. Let the elements 

aj,Q=p~''B{l+j^{p'^)) (modp"), j = 0,...,r-l 

he given. For each step k = 0, . . . ,u — 1 proceed as follows: 

The elements aj^k with j = 0, . . . ,rp — 1 have to be calculated successively by 

aj+r,k ^ i-iy^' E {-iraj+u,k (mod p"-*^). 

Set s/; = — ao,fc ^'[pi) (™od p) with < < p. The elements aj^^ which are divisible 
hy p are given by 

o^s^:+^^p,k = (modp), ^ = 0, ...,r-l. 

For k < u — 1 set 

aj,k+i = o.Sk+jp,k/p, i = 0, . . . , r - 1 

and go to the next step k-\-l, otherwise stop. Let (p, ti, . . . ,tn) S be the associated 
irregular pair with (p, I), then 

{p,ti, ...,tn,so,.. .,Su-i) e ^.;5^ 

is the only related irregular pair of order rn. 

Proof. Proposition 2.10 shows that 

j2(''){-ira,+,,0^0 (modp"). (5.1) 
i/=0 
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All elements of the sequence (aj.o)j>o can be calculated successively induced by the first 
elements aj^ with j = 0, . . . , r — 1. Using Corollary 2.5 with u = p'^(p[p'^), < k < u, 
then (5.1) becomes 



^ ('^V-ira^.+.p.^o ^ (modp"), (5.2) 



u=0 

whereby the sequence (aj+^pfe,o)M>o can also be calculated successively. Note that now 
j is fixed and /i runs. The sequences {ctj,k)j>o which we will consider are subsequences 
of (aj,o)j>o in a suitable manner. Essentially, these sequences are given by (5.2). The 
existence of these sequences and that they correspond to the related irregular pair of 
order n + k will be shown by induction on k for k = 0, . . . ,u — 1. Set In = I- By 
Proposition 2.7 and Theorem 3.1 there exist certain integers Sk and related irregular 
pairs of higher order for A; = 0, . . . , u — 1 where 

(p, In+k+l) G ^'n+fc+l> In+k+l = In+k + Sk (/^(p"'^'') with < Sk < P- 

Basis of induction k = 0: The sequence (aj.o)j>o is given by (5.1) and we have 

a,- = p-" B{ln + Mf')) (mod p"). 

Inductive step k k + 1: Assume true for k prove for k + 1. The elements aj^k with 
J = 0, . . . , r — 1 are given and the following elements are calculated by 

a,+r,k ^ (-1)'^+' E ( j (mod p^-") (5.3) 

i/=0 

up to index j = rp — 1. Proposition 2.7 provides 

Sk = -ao,fcA^|;^ (mod p) with < Sfc < p. (5.4) 

In the case k < u — 1, Lemma 2.6 additionally ensures that only Og^^jp^k = (mod p) 
for all j. Thus, we can define a new sequence by 

"j.fc+l = o;sk+jp,k/p 

^ p-{n+k+l) g^i^^^ ^ (^^ ^ ip)(^(p"+'')) 

^ p-("+'=+i)5(Wfc+i+jV(p"+'+')) (modp"-(^+i)) (5.5) 

for J = 0, ...,r — 1. By definition {paj^k+i)j=o,...,r-i is a subsequence of (aj,fe)j>o- 
Inductively (p'^^"^ aj^fc+i)j=o,...,r-i is a subsequence of (aj,o)j>o and satisfies (5.2) in 
a suitable manner and therefore also (5.3) considering case k + 1. On the other side. 
Congruence (5.5) shows that the new sequence also corresponds to the related irregular 
pair of order n + k + 1. 

Let {p,ti, . . . ,tn) € '^n'^ be the associated irregular pair with {p,l). Congruence 
(5.4) provides a unique integer for each step. Thus, {p,ti, . . . , tn, sq, ■ ■ ■ , Su-i) € "^1" 
is the only related irregular pair of order rn. □ 
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Remark 5.2 Unfortunately, the previous proposition has the restriction that for an 
irregular pair {p,l) € ^^'^ of order n and parameter r the following must hold that 

I > (r — l)n. 

Consider (691, 12) S ^'f. In this case one only could calculate related irregular pairs 
up to order 12. However, this restriction can be removed by shifting the index of the 
starting sequence (aj,o)j>o- Then shifting j ^ j + t yields 

l + tip{p"-) > (r - l)n 

enabling one to choose a greater value of r. In general, one has to proceed in the 
following way. Moving from a sequence (aj,fc)j>o to (aj,fc+i)j>o one has to determine 
elements aj k = (mod p). If one starts with a shifted sequence {a'j = {aj^k)j>t 
having elements aj± = (mod p) with < j < t, then the calculated sequence 
(cKjfc+i)j>o is also shifted in the index compared to {aj^k+i)j>o- This can be easily 
corrected by comparing the sequences resp. the resulting integers Sk with unshifted 
sequences calculated with a lower r' < r. In this case determining the integer Sk is 
better done by searching aj^k = (mod p) rather than calculating via (5.4). 

The main result can be stated as follows: irregular pairs of higher order can be 
determined with little effort by calculating a small number of divided Bernoulli numbers 
with small indices. We shall give another algorithm in terms of the p-adic zeta function 
(p^ I which produces a truncated p-adic expansion of X{p, i) ■ 

Proposition 5.3 Let {p,l) S with ^{p.i) 7^ 0. Let n be a positive integer. Initially 
calculate the values 

Cp,i,iik) = P'\p,iik) (modp"Zp) 

for k = 0, . . . ,n and 

A(p,i) = Cp,/, 1(0) - Cp,/, 1(1) (modpZp) 

where < '^(pj) < P- Set ti = 0. For each step r = 1, . . . , n proceed as follows: 
Calculate 

r 

^r = '^ Cp, 1, 1 (k) Tr^k (tr) (mod p^'Zp) 

k=0 

with the polynomials T^^k o,s defined in Theorem ^.^7. Then G p'''~^1,p. Set 

Sr+i = A^l^^ p^-'' (mod pZp) 

where < Sr+i < p- Set t^+i = tr + Sr+ip^~^ and go to the next step while r < n. 
Finally tn+i = X{p,i) (mod p"Zp). 

Proof. By Definition 4.20 we have Cp. «,i(^) = P^^Cpji^)- The value ^(pj) is given by 
its definition. Define t'^ = ipr-iiX{p,i)) for any r > 1 where the expansion of the zero 
of Cp,i is given by 

X{p,l) = S2 + S3P+ ... + Sr+lP"'^ + ... . 
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For now, let r E {1, . . . , n} be fixed. Theorem 4.27 provides 

r 

fc=o 

for s G Zp. From Definition 4.20 we liave 

Cp,Z,r(0) =P"''Cp,z(V'r--l(X(p,0)) =/~''Cp,;,i(4)- 
Proposition 4.22 shows that 

Cp,z,r(0) = A(p ;) (mod pZp). 

By construction we have t'^j^^ = t'^ + Sr+ip''^- Since t'^ = ti = we deduce by 
induction on r that t'^ = t^. and = Cp,l,i{t'r) (mod p''Zp) for 1 < r < n. This 
produces tn+i = i^n{X{p,i)) which equals to tn+i = X{p,i) (mod p^lip). Since 1 < r < n, 
we need the values Cp.l,i{k) (mod p^I^p) for A: = 0, . . . , n. □ 

Proposition 5.4 Xei n be an even positive integer, then 

5„ = (-l)t-l Jl ^r{p,n)+ordpn ^ "Q ^ 
p— Ifn p--l|n 

oo 

r(p, n) := J] #( n {(p, n mod (^(p^))} ). 
Here, as in Definition 2.1, x mod y denotes the least nonnegative residue of x modulo 

y- 

Proof. The trivial parts of the products above are given by (1.4), (1.5), and the 
sign. Thus, the product np-ifn^*^^^'"^ remains. From Definition 2.1 and Remark 2.2, 
the function r(p, n) follows by applying the maps X^, resp. Kummer congruences which 
results in a chain of related irregular pairs of descending order similar to (2.1). □ 

The previous proposition gives an unconditional representation of the Bernoulli 
numbers by means of the sets ■ Theorem 4.8 also gives a representation by zeros 
X{p,i) assuming the A-Conjecture. Of course, the problem of determining the occurring 
irregular prime factors remains open. On the other side, for instance, if one has cal- 
culated the first irregular pairs of order 10 for the first irregular primes pi, . . . ,pr like 
Table A. 3, then one can specify ad hoc all irregular prime powers pfj" with pi, < pr of 
Bn resp. C(l — n) up to index n = 4-10^^. Note that this lower bound is here determined 
by the first irregular prime 37 and order 10. 

Define for positive integers n and m the summation function of consecutive integer 
powers by Sn{m) = Y1^=q Many congruences concerning the function Sn are 
naturally related to the Bernoulli numbers. 
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Proposition 5.5 Let (p, I) G Assume that n = I (mod p — 1) where n > 0, then 

Bn _ Snip) 



n np 

Moreover 



(modp"). (5-6) 



^{p,i)=P [~rr[ ~) (™od^)) (5.7) 

with < A(p < p. 

Proof. Let n = I (mod p—1). The well-known formula of Sn, see [9, p. 234], is given 
by 

Snip) Bn^fn-1\ p2 " pk 

^ = - + V 1 )''''~'~3^^^[k-i)''-^'kikTT) 

where the equation is divided by n and p. Note that p > 37 and n > I > 12, because 37 
is the first irregular prime and -B12/I2 is the first divided Bernoulli number which has 
a numerator greater one. Now, the properties of (1.4) and (1.5) provide that Bn/n and 
Bn-2 are p-integral. For the other terms with Bn~k 7^ it follows that pBn-k is p- 
integral and ordp ip^~^ / ikik + 1))) > 2 by a standard counting argument. Therefore, 
Equation (5.8) is p-integral and holds (mod p^) whereas all terms of the right side 
vanish except for Bn/n. This gives Congruence (5.6). From Definition 2.3 we have 

pA(p,i) = B(Z+p-l)-^(/) (mod/). 

We can apply (5.6) in the congruence above. Reducing a^j-power and considering that 
/ + p- l = /-1^0 (mod p) finally yields (5.7). □ 

Looking at each line of Table A. 3, the product of the first three entries ^{p^i)-, si, 
and S2 are connected with the function Sn- Thus, one can easily verify these values. 

Proposition 5.6 Let ip,l) € ^'f with A(-p ;) 7^ 0. Let ip, 81,82) G be the related 
irregular pair of order two with I = si. Then 

A(p, 1)8182 = -p~^ Si ip) (mod p) . 
Proof. By Proposition 2.7 we have 

82 = -p^^yA^I^) (mod p). 
We then obtain by Proposition 5.5 that 

P82\pj) = -— = ^ (modp) 

which deduces the result since si = I < p. □ 
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Now, we shall give some reasons why a prediction or description of the occurrence 
of irregular prime factors of Bernoulli numbers seems to be impossible in general. For 
example, we have with an extremely small index n = 42 that 

1520097643918070802691 



observing that the numerator is a large irregular prime! As mentioned in Section 3, 
Bernoulli numbers Bn with index up to n = 152 have large irregular prime factors with 
30 up to 100 digits. This is even now the greatest mystery of the Bernoulli numbers! 

The connection with the Riemann zeta function C(s) via (1.1) leads to methods of 
calculating Bernoulli numbers directly in a fast and effective way, see [11, Section 2.7], 
noting that the main part of the calculation can be done using integers only. Let n be 
an even positive integer and \Bn\ = Un/Vn with {Un-,Vn) = 1 then (1.1) reads 

n} 1 — r 

^ ' p-l\n 

with Vn given by (1.4). Since C('^) ~^ 1 for ~^ "^n is a first approximation of the 
numerator [/„. Considering the decimal digit representation of Un and Tn, about n/3 
digits of the most significant decimal digits of C/„ and t„ are equal, see [11, Satz 2.7.9, 
p. 75] for a more precise statement and formula. Visiting 3^2 again, we observe 12 
identical digits: 

Ui2 = 1520 097 643 918070 802 691, 
r42 = 2 • 1806 • 42!/(27r)^2 ^ 1520 097 643 917 725172 488.7773. 

How can we interpret this result? A part of the most significant digits of the 
numerator of the Bernoulli number is determined in a certain way by all primes 
p < n + 1 and the reciprocal n-th power of vr. That the digits of vr are involved in the 
numerators of the Bernoulli numbers is quite remarkable. 



6 Connections with Iwasawa theory 

The A-Conjecture is directly connected with Iwasawa theory of cyclotomic fields over 
Q. Let p be an odd prime and /ipn be the set of p"-th roots of unity where n is a positive 
integer. For the cyclotomic field Q(/ipi) let Q(/ip")''' be its maximal real subfield. The 
class number hp = h{Q{iJ.p)) can be factored by hp = hp h^ where h^ = h{Q{ij,p)~^) 
and hp is the relative class number introduced by Kummer. Define 

Bi^ru = - auj'^ia) 

a=l 

as the generalized Bernoulli number assigned to the Teichmiiller character uj. This 
character is defined by w : Z* ^ Z* and uj{a) = a (mod p'^p) for a S Z* giving the 
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(p — l)-th roots of unity in Qp. We have for even m > and p — \ \ m the following 
relation, see [23, Corollary 5.15, p. 61]: 

-^i,!^™-! = B{m) (mod pZp). 

For the detailed theory, especially of Iwasawa invariants and cyclotomic Zp-exten- 
sions, see Washington [23] and Greenberg [8]. The results of Iwasawa, Ferrero and 
Washington, Vandiver and Kummer provide the following theorem, see [23, Cor. 10.17, 
p. 202]. 

Theorem 6.1 Let p he an irregular prime. Assume the following conditions for all 
irregular pairs (p, I): 

(1) The conjecture of Kummer- Vandiver holds: p\hp, 

(2) The Kummer congruence does not hold (mod p^): B{1 +p—l) ^ B{1) (mod p^), 

(3) The generalized Bernoulli number is not divisible by p^ : Bi^i-i ^ (mod p^ljp). 
If these are satisfied, then 

ordp /i(Q(/Xpn)) = i{p)n for all n > 1. 

All conditions of the theorem above hold for all irregular primes p < 12 000 000 as 
verified in [2]. In the case of a regular prime p the formula of the theorem above is also 
valid, because then we have i{p) = and 

p \ hp = h{Q{pp)) <^ P^h{Q{fip^)) for all n>l 

as follows, e.g., from Iwasawa theory. To get another point of view we can exchange 
two conditions of the previous theorem by our results. Conditions equivalent to those 
of Theorem 6.1 are as follows: 

(2') The A-Conjecture holds: A(p^;) ^ 0, 

(3') A special irregular pair of order two does not exist: (p, 1,1 — 1) ^ ^'2"- 

Now, the A-Conjecture with its consequences gives a significant reason to believe 
that the condition (2) resp. (2') may hold in general. We still have to show that the 
condition (3') is equivalent to the condition (3). 

Proposition 6.2 Let {p,l) G "^^^ . Then 

Bi,uji-i = -Bi+(p„i)(i_i) (mod p^Zp) 

and 

= (mod p%,) ^ {p, 1,1 -l)e $'2"- 

To prove this result, we first need some properties of the Teichmiiller character 
uj. Since u!{a) is defined by uj{a) = lim„_too o''" in the following lemma is easily 
derived. 
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Lemma 6.3 Let a, p be integers with p an odd prime and < a < p. Then 

uj(a) = (mod p^Zp) and uj{a) = aF+p{a^ — a) (mod p'^Zp). 

Proof of Proposition 6.2. Since a;'~^(a) = Lij(a'~^), we have 

p— 1 p— 1 

a=l a=l 

Using Lemma 6.3 we obtain 

aw(a'~i) EE a^^('-i)+i + p (aP('~i)+i - a') (mod p=^Zp). 

From the definition of 5„ and Proposition 5.5 we have 

P^iy-i = 5'p(;_i)+i(p) +p5p(j_i)+i(p) -pSiip) = pBi^^p_-^)^i^]^) (mod /Zp). 

Since p(/ - 1) + 1 = / + (p - 1)(Z - 1) = / (mod p - I) and (p, Z) G ^'i", only the first 
term S'p(;_i)+i(p) does not vanish. Note that p\ p{l — 1) + 1, therefore we get 

= i3,+(p_i)(,_i) = B{1 +{p- - 1)) (mod p") 

if and only if {p, 1,1-1) e ^'f. □ 

Remark 6.4 From Propositions 5.5 and 5.6, the conditions '^(pj) 7^ and {p, I, l — l) ^ 
^'2'^ are equivalent to the system 

I Si+p.i{p) - (/ - 1) Slip) ^ (mod p^), 
I Si+p.i{p) -{1-2) Slip) # (mod p^). 

7 The singular case 

In Section 4 we have derived most of the results assuming the A-Conjecture. Theorem 
4.8 conjecturally describes a closed formula for ^(1 — n) by zeros X(p,i)- The following 
theorem gives an equivalent formulation for the Bernoulli numbers. 

Theorem 7.1 Let n be an even positive integer. Under the assumption of the A- 
Conjecture we have 

Bn=i-i)^-' n \< n \pixip,i)-^,)\p' 

p-lfn (p,z)e*i^^ p-l\n 

l = n (mod p— 1) 
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Proof. We have to modify the formula of Theorem 4.8. The product formula gives 
1= |n|p = jnjoo \n\p |?7,|p . (7.1) 

pGPU{oo} p-l|n p-l\n 

Since —Bn/n = C,{\ — n), the proposed formula follows easily. □ 

To get an unconditional formula for Bn resp. C(l — n) we have to include the case 
of a singular A(p However, no such singular ^(pj) has been found yet. Theorem 
3.2 describes the more complicated behavior of related irregular pairs of higher order 
in the singular case which can be described by a rooted p-ary tree, see Diagram 3.3. 

Let {p,l) £ ^i'^ with a singular A(p /). We construct the rooted p-ary tree of 
related irregular pairs of higher order which is a consequence of Theorem 3.2. Each 
node contains one related irregular pair of higher order. Note that these pairs are not 
necessarily distinct. We denote this tree by T^^ assigned to the root node (p, I). 

The tree T^^ has the property that each node of height r lies in A tree 

-^(p ~ {(^''Ol is called a trivial tree having height 0. If the tree T^'^ is of height 
> 1, then it contains the root node {p,l) and its p child nodes {p,l + j^p{p)) € \I'2" for 
j = 0,...,p-l. 

In the nonsingular case, we have a zero of the p-adic zeta function. In contrast 
the singular case does not guarantee that related irregular pairs of higher order exist 
at all. The discovery of a singular A(p ;) is not incompatible with Theorem 7.1 but 
the formula becomes more complicated, because we then have to consider the complete 
tree T^P . By combining both cases we obtain an unconditional formula which is given 
by the following theorem. Recall Definition 2.1. 

Theorem 7.2 Let n be an even positive integer, then 

Bn = (-1)^-^ n i^ip' n \pix(p,i)-^i)\;' 

p-l\n {p,OG*i", A(p_;j#o 

Z = n (mod p—1) 

(p,^)e*f^A(p_;)=o p-l|n 

l = n (mod p — 1) 

with the height h'^^ of n defined by 

h^j^in) = max{height((p,r)) : {pJ') G T^,i) n {(p,n mod (^(p'^))}.>i} . 
Moreover, h^^^ /)(^) = tree T^^ is trivial. 

Proof. The case / is already covered by Theorem 7.1. Next we assume that 

A(p,z) = with a given tree where n = / (mod p — 1). As a consequence of the 

construction of T^^ and Remark 2.2, we have to determine the maximal height of a 
node {p, lyj) G T^^ fl ^''J'' where {p, l^j) = {p, n mod Lp{jp^)). The root node (p, I) has 
height 0, so the exponent equals 1 + /i9 An). 
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If the tree has the height > 1, then (p, n mod ^{p^)) G T'^ ^y, this imphes that 

^(p O^"'-' — 1- ^ trivial tree T^P imphes that h^^ ^^(n) = 0. Conversely, if h^^ = Oi 
then the height of the tree T^^ must be zero, otherwise we would get a contradiction. 

Corollary 7.3 Let n he an even positive integer, then 

c(i-n)=(-i)t n iMx(„o-^)i;^ n p n ^ 

(p.i)e*i",A(p_,)7^o {p,Oe*i",A(p ,)=o p-l|n 

/ =n (mod p — 1) l = n {mod p—1) 

with as defined above. 

Proof. This is a reformulation of Theorem 7.2 by ("(1 — n) = —Bn/n and (7.1). □ 



8 An extension of Adams' theorem 

Let n be an even positive integer. The trivial factor of Bn, given by (1.5), is a conse- 
quence of the implication, known as Adams' theorem, that 

p^ I n with p — l\n p' \ Bn 

where p is a prime and r is some positive integer. It was, however, never proved by 
Adams. In 1878 he computed a table of Bernoulli numbers i?2m for m < 62. On the 
basis of this table he conjectured^ that p j n implies p \ Bn for primes with p— 1 f n, see 
[1]. Note that the property that Bn/n is a j5-integer for p — 1 f n is needed to formulate 
the Kummer congruences (1.3). The case r = 1 of these congruences was proved by 
Kummer [15] earlier in 1851. 

By Theorem 7.2 and the definitions of h^^ and X{p,i)^ we can state an extended 
version of Adams' theorem. We introduce the following notation. We write p'' \\ n 
when p^ I n but p^~^^ f n, i.e., r = ordp n. 

Theorem 8.1 Let n be an even positive integer. Let p be a prime with p^ \ \ n, r > 1, 
and p — 1 \ n. Let I = n (mod p—1) with < / < p — 1. Then p^'^^ \\ Bn with the 
following cases: 

(1) If p is regular, then (5 = 0, 

(2) If p is irregular with {p,l) ^ "^i^, then 5 = 0, 

(3) Ifp is irregular with {p, I) £ ^'f , A(p ;) / 0, then 6=1 + ordp (x(p,;) - fEj), 

(4) If p is irregular with {p, I) € ^i^, '^(p.i) = 0, then 6=1 + h^^ ;)("-)• 
Additionally, in case (3) resp. (4), if {p,l,l) ^ ^2^^, then 5=1, otherwise 6 >2. 



t "... I have also observed that if p be a prime factor of n which is not likewise a factor of the denominator 
of the nth number of Bernoulli, then the numerator of that number will be divisible by p. I have not 
succeeded, however, in obtaining a general proof of this proposition, though 1 have no doubt of its truth." 
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Proof. We have to consider the formula of Theorem 7.2. The first product yields 
I Bn- Only the second resp. third product can give additional p- factors. Therefore 
case (1) and (2) are given by definition. We can now assume that {p, I) € ^''i^''. 

Case (3): A nonsingular ^(pj) provides 

6 = ordp \p{x{p,l) - ^ = 1 + ordp {x{p,l) - 

By assumption n = p'^'n' with some integer n'. We have to evaluate 

d = ordp (x(p,i) - ^) = ovdp{px{p,i) - X(p,z) + l-p''n). 
Since r > 1, we p-adically obtain 

(p, I, I) G ^ X{p,l) =l + S3p+... ^ d>l. 

Conversely, {p, I, I) ^ "if 2^ yields d = 0. 

Case (4): A singular A(p provides 6 = 1 + hP^^ The definition of and 

Theorem 7.2 show that 

(p,/,/) ^ $'2" the tree T(° is trivial ^('p,o(") = ^■ 

Conversely, (p, /,/) S "^'f^ yields /i^^ ^^(n) > 1. □ 

So far, no (p, Z, /) E ^'2'^ has been found. The following corollary theoretically shows 
examples where 5 is arbitrary large. 

Corollary 8.2 Assume that {p,l, . . . ,1) G ^r+i exists with some r > 1. Set n = Ip^ . 

Then we have p^ \ \n and p^^+i j j^^^ j_g_^ ^ > _|_ x. 

Proof. By Definition 2.11 {p,n) G ^I/Jf^^ is associated with {p,l,. .. ,1) G ^'J.+i, since 
n = /p'' = J2lt\ Mp""^)- Thus, p''+^ I Bn/n and finally p^^'+i \ B^. □ 

Remark 8.3 As mentioned above, Johnson [10] calculated the now called irregular 
pairs (p, si,S2) G ^2" of order two for p < 8000. He also proved that {p,l,l) ^ ^2" 
resp. si 7^ S2 in that range. In a similar manner, the nonexistence of irregular pairs 
(p, Z, Z — 1) of order two plays an important role in Iwasawa theory as seen in Section 6. 
One may conjecture that no such special irregular pairs (p, I, I) and (p, 1,1 — 1) of order 
two exist. But there is still a long way to prove such results, even to understand 
properly which role the zeros X{p, i) play. Now, we have the relation 

(p,Z,0^$r ^ P^\B{lp) ^ p^\Bip. 

Yamaguchi [24] also verified by calculation that p^ f Bip for all irregular pairs (p, I) 
with p < 5500, noting that this was conjectured earlier by Morishima in general. The 
condition p^ | Bip is related to the second case of FLT, see [23, Thm. 9.4, p. 174]. 
See also [23, Cor. 8.23, p. 162] for a different context. Under the assumption of the 
conjecture of Kummer-Vandiver and that no (p, /,/) G ^'2'^ exists, the second case of 
FLT is true for the exponent p. For details we refer to the references cited above. 
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The converse of Adams' theorem does not hold, but one can state a somewhat differ- 
ent result which deals with the common prime factors of numerators and denominators 
of Bernoulli numbers with indices close to each other. 

Proposition 8.4 Let S = {2, 4, 6, 8, 10, 14} be the set of all even indices m where the 
numerator of \Bm/m\ equals 1. Write Bn = An/Vn with (yl^, Vn) = 1. Let k, n he even 
positive integers with k ^ S and n — k > 2. Then 

D = {An,Vn-k) implies D \ n. 

Moreover, if D > 1 then D = Pi ■ ■ ■ Pr with some r > 1. The primes pi, . . . ,pj. are 
pairwise different and py \ V^, Pu \ Bnjn for i' = 1, . . . ,r. 

Proof. Assume that D > 1. We then have D = pi ■ ■ - pr with some r > 1, since 
Vn-k is squarefree by (1.4). Let i/ € {1, . . . , r}. Since pi, \ An and pi, \ Vn-k, we have 
Pi, — l]n and Pu — 1 \ n — k. From this we can deduce that p^ — 1 \ k and consequently 
that Pu \Vk- Next we assume that p^ \ n or p^ \ Bn/n. Note that p^ \ An and Pu \ n 
imply that p^ \ Bn/n, but not conversely. We can use the Kummer congruences (1.3) 
to obtain that 

Bn -B/j 

0= — = — (modp^), 
n k 

since n = k (mod p^ — 1). By the definition of the set S we have 

(modpu) (8.1) 

which yields a contradiction. This shows that p^ \ n and p^ \ Bn/n. Finally it follows 
that D \n. □ 

Now, the set S cannot be enlarged, because (8.1) does not hold in general for 
numerators having prime factors. For example, let p = 691 and n = 12 + (p — 1) = 702, 
then we have p \ -B12/I2 and D = {An, V^-12) = pc\n with some c > 1. On the other 
hand, one trivially obtains for k G S, p prime with p — 1 \ k, and n = kp infinitely 
many examples of > 1. In the following proposition. Proposition 8.4 plays a crucial 
role. Recall the definition of S'„(m). 

Proposition 8.5 Let n,m be positive integers with even n. For r = 1, 2 we have 

m^^^ I Sn{rn) <J=^ rnJ' \ Bn- 



Proof. We can assume that m > 1, since m 
B2 = and that f |m(m — l)(2m — 1) = 
that n > 4. We have, see (5.8), that 



= 1 is trivial. The case n = 2 follows by 
= 5*2 (?n,) for m > 1. For now we assume 
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By (1.4) and the cases Bq = I and Bi = the denominator of all nonzero Bernoulli 
numbers is squarefree. For each prime power factor \\ m and k where -Bn-fc 7^ 
{2 < k < n) we have 

L j^n-fc-^^j >s(fc + l)-l-ordp(A; + l) > As (8.3) 

with the cases: (1) A = 1 for A; > 2, p > 2, (2) A = 2 for A; > 2, p > 5, and (3) A = 3 
for /c > 4, p > 5. The critical cases to consider are p = 2,3,5 and s = 1. Now, we are 
ready to evaluate (8.2) (mod m'") for r = 1, 2. Write i?„ = An/Vn with {An, Vn) = 1. 

Case r = 1: Assume that {m, Vn) > 1. By (8.3) (case A; > 2, p > 2) we obtain 

Sn{m) = Bn m = — ^ m ^ (mod m). 

Therefore, (m, Vn) = 1 must hold which implies 2 \ m, 2> \ m, and p > 5. Hence, by 
(8.3) (case A; > 2, p > 5), we can write Sn{in) = Bn m (mod m?). This yields 

I Sn{m) <^ m\Bn. (8.4) 

Case r = 2: We have m \ Bn and (m, 6) = 1, because either | Bn or rrfi \ Sn{fn) 
is assumed. The latter case implies | Sn{m) and therefore m \ Bn by (8.4). Since 
|yl4| = 1, we can assume that n > 6. We then have Bns = and we can apply (8.3) 
(case A: > 4, p > 5) to obtain 

Sn{m) =Bnm+ '^^"^ ~ ^^^""-^ (mod jn^). (8.5) 

6K-2 

Our goal is to show that the second term of the right side of (8.5) vanishes, but 
the denominator Vn-2 could possibly remove prime factors from m. Proposition 8.4 
asserts that (An, Vn-2) \ n. We also have {m,Vn~2) \ n since m \ Bn- This means that 
the factor n contains those primes which Vn-2 possibly removes from m. Therefore 
the second term of (8.5) vanishes (mod m^). The rest follows by ^^(m) = Bnm = 
(mod m^). □ 

One cannot improve the value r in general. Choose p = 37 and / = 37580. Since 
{p,l) G we have p^ \ Bi, but p^ \ Si{p) which was checked with Mathematica. 

Example 8.6 

(1) We have B42 = 1520097643918070802691/1806. Since the numerator yl42 is a 
large irregular prime, we obtain for m > 1 that 

I 5'42(?Ti) m = 1520097643918070802691. 

(2) We have A50 = 5^ • 417202699 • 47464429777438199 and V"48 = 2 • 3 • 5 • 7 • 13 • 17. 
Hence, for m > 1 we have 

I 55o(m) <;=^ m = 5. 
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A Calculations 

Table A.l Bn and B„/n. 



n 





1 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 




1 


1 

2 


1 
6 


1 

30 


1 

42 


1 

30 


5 

66 


691 
2730 


7 
6 


3617 
510 


43867 
798 


174611 
330 


s„ 






1 


1 


1 


1 


1 


691 


1 


3617 


43867 


174611 


n 






12 


120 


252 


240 


132 


32760 


12 


8160 


14364 


6600 



Table A. 2 Calculated irregular pairs of order 100 of primes 37, 59, and 67. 



Case p = 37. 



Zeros of the sequence (sv) occur at index 19 and 81. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 





32 


7 


28 


21 


30 


4 


17 


26 


13 


32 


10 


35 


27 


36 


32 


10 


21 


9 


11 





1 


20 


13 


6 


8 


10 


11 


10 


11 


32 


13 


30 


30 


10 


6 


8 


2 


12 


1 


8 


2 


5 


3 


40 


10 


19 


8 


4 


7 


19 


27 


33 


29 


29 


50 


11 


2 


23 


8 


34 


5 


8 


35 


35 


13 


60 


31 


29 


6 


7 


22 


13 


29 


7 


15 


22 


70 


20 


19 


29 


2 


14 


2 


2 


31 


11 


4 


80 





27 


8 


10 


23 


17 


35 


15 


32 


22 


90 


14 


7 


18 


8 


3 


27 


35 


33 


31 


6 




Case p = 59. 

Zeros of the sequence (s^) occur at index 31 and 95. 


s„ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 





44 


15 


25 


40 


36 


18 


11 


17 


28 


58 


10 


9 


51 


13 


25 


41 


44 


17 


43 


35 


21 


20 


10 


21 


38 


9 


12 


40 


43 


45 


30 


41 


30 





3 


25 


34 


49 


45 


9 


19 


48 


57 


40 


11 


13 


29 


28 


44 


41 


37 


33 


29 


43 


50 


8 


57 


12 


48 


15 


15 


53 


57 


16 


51 


60 


16 


54 


30 


9 


26 


8 


49 


22 


58 


11 


70 


42 


28 


36 


33 


45 


24 


32 


18 


12 


29 


80 


45 


40 


27 


19 


40 


41 


11 


42 


49 


35 


90 


41 


57 


54 


33 





34 


34 


49 


6 


31 




Case p = 67. 

Zeros of the sequence (sv) occur at index 23 and 85. 


s„ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 





58 


49 


34 


42 


42 


39 


3 


62 


57 


19 


10 


62 


10 


36 


14 


53 


57 


16 


60 


22 


41 


20 


21 


25 





56 


21 


24 


52 


33 


28 


51 


30 


34 


60 


8 


47 


39 


42 


33 


14 


66 


50 


40 


48 


45 


28 


61 


50 


27 


8 


30 


59 


32 


50 


15 


3 


1 


54 


12 


30 


20 


14 


12 


10 


60 


49 


33 


49 


54 


13 


26 


42 


8 


58 


12 


70 


63 


19 


16 


48 


15 


2 


13 


1 


23 


2 


80 


44 


64 


25 


40 





16 


58 


44 


31 


62 


90 


47 


61 


46 


9 


2 


50 


1 


62 


34 


31 



38 



Table A. 3 Calculated irregular pairs of order 10 of primes below 1000. 



(p,0 


^(P,0 


Sl 


S2 


S3 


S4 


S5 


S6 


S7 


sg 


s<j 


sio 


(37,32) 


21 


32 


7 


28 


21 


30 


4 


17 


26 


13 


32 


(59,44) 


26 


44 


15 


25 


40 


36 


18 


11 


17 


28 


58 


(67,58) 


21 


58 


49 


34 


42 


42 


39 


3 


62 


57 


19 


(101,68) 


42 


68 


57 


57 


45 


60 


16 


10 


47 


53 


88 


(103,24) 


54 


24 


2 


87 


55 


47 


3 


72 


4 


45 


52 


(131,22) 


25 


22 


93 


26 


43 


74 


109 


80 


5 


55 


14 


(149,130) 


79 


130 


74 


68 


10 


94 


16 


122 


70 


110 


10 


(157,62) 


48 


62 


40 


145 


67 


29 


69 





87 


89 


21 


(157,110) 


51 


110 


73 


3 


58 


9 


114 


118 


21 


1 


11 


(233,84) 


132 


84 


173 


164 


135 


146 


127 


10 


36 


108 


230 


(257,164) 


188 


164 


135 


174 


30 


203 


161 


193 


142 


68 


126 


(263,100) 


87 


100 


198 


139 


151 


106 


202 


99 


202 


251 


163 


(271,84) 


179 


84 


5 


14 


239 


8 


233 


43 


28 


57 


170 


(283,20) 


15 


20 


265 


115 


171 


137 


251 


118 


132 


246 


265 


(293,156) 


93 


156 


230 


75 


289 


47 


247 


98 


100 


141 


27 


(307,88) 


205 


88 


70 


234 


51 


173 


104 


140 


140 


107 


201 


(311,292) 


277 


292 


204 


183 


9 


260 


183 


214 


254 


2 


151 


(347,280) 


106 


280 


113 


250 


150 


307 


264 


145 


177 


101 


156 


(353,186) 


301 


186 


190 


147 


13 


34 


171 


106 


304 


190 


102 


(353,300) 


161 


300 


181 


300 


314 


327 


67 


26 


113 


18 


336 


(379,100) 


276 


100 


242 


277 


88 


236 


225 


22 


221 


54 


26 


(379,174) 


82 


174 


364 


216 


20 


128 


277 


134 


257 


164 


31 


(389,200) 


48 


200 


354 


33 


371 


189 


29 


219 


44 


11 


319 


(401,382) 


376 


382 


263 


126 


213 


197 


170 


320 


107 


297 


331 


(409,126) 


180 


126 


389 


343 


247 


322 


24 


187 


75 


91 


179 


(421,240) 


396 


240 


351 


141 


36 


169 


124 


164 


342 


365 


156 


(433,366) 


284 


366 


406 


342 


372 


234 


21 


328 


346 


279 


155 


(461,196) 


281 


196 


423 


121 


233 


61 


353 


421 


414 


350 


92 


(463,130) 


78 


130 


376 


404 


124 


420 


63 


438 


185 


124 


18 


(467,94) 


118 


94 


219 


393 


264 


70 


75 


254 


361 


332 


157 


(467,194) 


269 


194 


283 


329 


154 


419 


170 


152 


78 


304 


326 


(491,292) 


456 


292 


218 


299 


225 


362 


461 


37 


65 


203 


228 


(491,336) 


103 


336 


260 


15 


41 


381 


66 


376 


391 


209 


305 


(491,338) 


475 


338 


59 


160 


106 


105 


33 


346 


158 


314 


233 


(523,400) 


497 


400 


36 


230 


180 


431 


235 


114 


104 


152 


399 


(541,86) 


211 


86 


436 


29 


482 


424 


74 


212 


259 


419 


287 


(547,270) 


348 


270 


458 


536 


35 


521 


413 


88 


545 


44 


537 


(547,486) 


139 


486 


100 


4 


33 


153 


282 


467 


233 


482 


17 


(557,222) 


153 


222 


549 


505 


399 


472 


49 


20 


81 


279 


513 


(577,52) 


452 


52 


309 


416 


274 


56 


20 


476 


164 


309 


19 


(587,90) 


286 


90 


109 


344 


244 


53 


93 


454 


292 


291 


547 


(587,92) 


319 


92 


213 


332 


470 


36 


479 


508 


134 


323 


275 


(593,22) 


331 


22 


188 


388 


541 


576 


371 


26 


586 


40 


514 


(607,592) 


435 


592 


369 


428 


162 


503 


358 


484 


411 


67 


267 


(613,522) 


57 


522 


549 


451 


318 


312 


243 


38 


265 


552 


215 


(617,20) 


289 


20 


384 


107 


161 


281 


358 


64 


604 


336 


326 


(617,174) 


317 


174 


546 


83 


114 


484 


121 


229 


335 


597 


570 


(617,338) 


312 


338 


419 


570 


496 


63 


247 


46 


604 


464 


134 


(619,428) 


121 


428 


457 


363 


526 


36 


179 


79 


170 


485 


47 


(631,80) 


139 


80 


146 


468 


175 


34 


249 


169 


26 


498 


528 


(631,226) 


221 


226 


338 


510 


318 


581 


572 


363 


422 


111 


405 


(647,236) 


318 


236 


480 


525 


205 


103 


205 


620 


394 


553 


25 


(647,242) 


94 


242 


487 


519 


49 


109 


373 


451 


586 


250 


57 


(647,554) 


209 


554 


558 


568 


174 


579 


545 


5 


377 


242 


81 


(653,48) 


363 


48 


154 


558 


439 


300 


59 


541 


242 


205 


47 


(659,224) 


200 


224 


140 


131 


396 


158 


367 


79 


256 


620 


615 


(673,408) 


325 


408 


26 


64 


257 


158 


213 


430 


659 


144 


600 


(673,502) 


585 


502 


293 


198 


436 


506 


441 


27 


89 


416 


407 



39 



ipj) 




Sl 


S2 


S3 


Si 


S5 


S6 


ST 


sg 


S9 


sio 


(677,628) 


440 


628 


504 


457 


324 


461 


88 


532 


653 


89 


244 


(683,32) 


477 


32 


266 


20 


625 


119 


190 


13 


190 


222 


214 


(691,12) 


611 


12 


496 


104 


197 


607 


590 


303 


96 


461 


152 


(691,200) 


592 


200 


496 


333 


578 


93 


160 


436 


611 


215 


278 


(727,378) 


398 


378 


683 


722 


169 


391 


150 


694 


210 


228 


130 


(751,290) 


164 


290 


481 


37 


181 


27 


31 


71 


8 


36 


164 


(757,514) 


554 


514 


364 


164 


375 


7 


720 


750 


273 


592 


643 


(761,260) 


462 


260 


729 


680 


274 


188 


464 


183 


283 


52 


235 


(773,732) 


517 


732 


147 


306 


278 


370 


412 


89 


340 


637 


223 


(797,220) 


375 


220 


369 


279 


501 


300 


168 


530 


534 


747 


268 


(809,330) 


88 


330 


52 


743 


100 


336 


157 


759 


348 


43 


736 


(809,628) 


18 


628 


773 


629 


623 


160 


494 


339 


244 


463 


274 


(811,544) 


381 


544 


424 


100 


346 


749 


624 


220 


410 


313 


62 


(821,744) 


704 


744 


621 


319 


498 


427 


50 


21 


237 


305 


809 


(827,102) 


105 


102 


164 


443 


469 


568 


671 


183 


372 


512 


464 


(839,66) 


269 


66 


135 


305 


36 


40 


659 


431 


326 


591 


293 


(877,868) 


480 


868 


554 


279 


714 


821 


520 


76 


565 


104 


22 


(881,162) 


789 


162 


372 


330 


89 


244 


27 


229 


418 


438 


89 


(887,418) 


611 


418 


76 


698 


835 


872 


130 


319 


217 


439 


573 


(929,520) 


607 


520 


433 


27 


711 


366 


902 


838 


7 


351 


805 


(929,820) 


706 


820 


749 


156 


59 


913 


480 


432 


114 


129 


491 


(953,156) 


24 


156 


720 


516 


620 


229 


251 


77 


805 


689 


477 


(971,166) 


715 


166 


538 


594 


897 


509 


355 


749 


180 


174 


96 



Table A. 4 Calculation: n = 1, p = 37, I = 32, {p, I) G 



j 


Index 


a j (mod p"^) 


= (mod p'^) 


(mod p-^) 


Aaj (mod p^) 





32 


3941/2720 


42144 


45827 


650 


1 


68 


2587/15 


37318 


49934 


650 


2 


104 


3821/1272 


36599 


30768 


650 


3 


140 


6497/7198 


16714 





Using Proposition 2.10 with r = 3 and (r — l)n = 2 yields s = 1043 (mod p^) 
and = 32 + sip{p) = 37580. We obtain (37,284) G ^-^^ (37,37580) e ^f, and 
(37,32,7,28) G 

Table A.5 Calculation: n = 3, p = 37, / = 37580, {p, I) G "^'f. 



j 


Index 


Oj (mod p'^) 


= (mod p'^) 


Aa, (mod p^) 





37580 


11241/22913 


24645 


45827 


1 


86864 


49609/46188 


19819 


45827 


2 


136148 


5261/24 


14993 





Using Proposition 2.9 yields s = 6607 (mod p^) and Iq = 37580 + sip{p^) = 

325656968. We obtain (37,325656968) G (37,55777784) G (37,1072544) G 
, and (37, 32, 7, 28, 21, 30, 4) G 

Table A. 6 Calculation: n = 3, p = 37, / = 37580, {p, I) G ^jj"'. 



j 


Index 


Qj (mod p^) 


= (mod p^) 





37580 


3791602112159/3307480 


45520991695194 


1 


86864 


1046892158059/484258896735 


47985230204445 


2 


136148 


13280633201029/15 


70198303437443 


3 


185432 


8822143378793/98280020 


73479320052104 



Using Proposition 5.1 with r = 4 and (r — l)n = 9 yields the sequence 
21, 30, 4, ... , 27 which provides (37, 32, 7, 28, 21, 30, 4, 17, 26, 13, 32, 35, 27) G ^g. 
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Note that Tables A. 2 and A. 3 were calculated with smallest possible indices of 
the Bernoulli numbers using Proposition 5.1; they agree with these results above. 
Additionally, the results were checked by Corollary 4.23 and Proposition 5.3. The 
program calcbn [11, Section 2.7] was used to calculate these large Bernoulli numbers 
extremely quickly. 
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